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Abstract 

This  paper  describes  a  comparative  investigation  of  friction-compensating  con¬ 
trol  strategies  designed  to  improve  low-velocity  position  tracking  performance 
for  servomechanisms.  Several  control  methods  are  considered  including  adap¬ 
tive  control  and  estimation-based  control.  Additionally,  the  various  controller 
designs  incorporate  different  friction  models  ranging  from  classical  friction  and 
Stribeck  friction  to  the  less  popular  Dahl  friction  model.  This  investigation  of 
friction  models  is  motivated  by  the  fact  that  there  is  little  consensus  in  the  lit¬ 
erature  on  how  best  to  model  friction  for  dynamic  friction  compensation.  The 
control  strategies  are  compared  in  an  extensive  test  program  involving  posi¬ 
tion  tracking  experiments  on  a  direct-drive  dc  motor.  This  effort  addresses  the 
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current  lack  of  comparative  experimental  results  on  friction  compensation.  The 
results  show  that  the  adaptive  and  estimation-based  controllers  outperform  more 
traditional  linear  controllers.  The  experiments  also  yield  insight  into  the  appro¬ 
priateness  of  the  different  friction  models  under  the  tested  operating  conditions. 
In  particular,  the  Dahl  model  is  observed  to  provide  a  reliable  representation  of 
friction  behavior  near  zero  velocity. 


1  Introduction 

Recent  growth  in  the  number  and  variety  of  robotics  applications  has  led  to  a 
demand  for  increased  precision  in  robotic  manipulation.  However,  robotic  ma¬ 
nipulators  must  contend  with  friction  which  poses  a  serious  challenge  to  precise 
manipulator  control.  Specifically,  failing  to  compensate  for  friction  can  lead  to 
tracking  errors  when  velocity  reversals  are  demanded  and  oscillations  when  very 
small  motions  are  required.  To  compensate  for  friction  it  is  best  to  have  some 
knowledge  of  the  structure  of  friction,  yet  there  is  little  overall  agreement  in 
the  literature  on  how  best  to  model  friction.  Further,  friction  compensation  is 
complicated  by  the  fact  that  friction  parameters  vary  with  temperature  and  age. 

Traditionally,  control  engineers  have  used  open-loop  smoothing  techniques, 
such  as  dither  and  pulse-width  modulation,  to  compensate  for  friction  in  me¬ 
chanical  systems.  However,  these  techniques  have  disadvantages,  for  example, 
dither  can  cause  mechanical  problems  such  as  fatigue  by  exciting  vibrations  in 
manipulators. 

As  an  alternative  to  these  techniques,  researchers  in  control  engineering 
have  recently  considered  adaptive  and  estimation-based  control  techniques  for 
compensation  of  friction  in  mechanical  systems  [Gilbart  and  Winston,  1974, 
Walrath,  1984,  Craig,  1988,  Canudas  et  al.,  1986].  Among  the  compensators 


proposed  in  the  literature,  a  variety  of  friction  models  are  assumed.  Much  of 
this  work  shows  experimental  evidence  that  a  particular  friction  model  together 
with  a  suitable  compensation  technique  improves  system  performance.  How¬ 
ever,  since  each  research  team  performed  different  experiments,  there  is  no  easy 
way  to  compare  the  relative  effectiveness  of  the  different  control-technique  and 
friction-model  combinations. 

This  paper  provides  experimental  results  that  indicate  the  relative  effec¬ 
tiveness  of  five  different  friction-compensating  controllers.  Of  the  five  con¬ 
trollers  tested,  two  are  modified  versions  of  adaptive  controllers  designed  by 
[Gilbart  and  Winston,  1974]  and  [Craig,  1988],  respectively.  Each  assumes  the 
classical  model  of  friction.  However ,  in  the  present  study  both  adaptive  con¬ 
trollers  have  been  upgraded  to  include  more  detailed  friction  elements  such  as 
asymmetries  and  Stribeck  friction. 

The  third  controller  tested  is  a  modified  version  of  the  estimation-based 
friction-compensating  controller  of  [Walrath,  1984].  This  controller  design  is 
based  on  the  Dahl  friction  model  which  predicts  a  first-order  dynamic  model 
of  friction  as  a  function  of  displacement  with  a  time  constant  that  is  a  linear 
function  of  velocity.  Whereas  Walrath  was  unable  to  experimentally  derive  this 
linear  time  constant  function,  the  experiments  described  in  this  paper  success¬ 
fully  verify  Dahl’s  prediction.  Additionally,  an  original  stability  proof  for  the 
estimation-based  controller  that  uses  the  passivity  formalism  is  outlined  in  this 
paper. 

The  fourth  controller  tested  is  a  linear  controller  with  dither.  This  controller 
is  included  as  a  sample  smoothing  controller.  The  benchmark  for  the  test  pro¬ 
gram  is  a  conventional  linear  controller  with  optimized  proportional,  integral, 
and  derivative  (PID)  gains. 
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The  experimental  program,  designed  for  the  hardware  available  in  the  lab¬ 
oratory,  provides  a  realistic  servomechanism  control  problem.  The  subject  of 
the  program  is  a  direct-drive  brush-type  dc  motor,  digitally  controlled  by  means 
of  an  IBM  AT  personal  computer  (PC).  The  same  series  of  position  tracking 
experiments  are  performed  on  the  motor  with  each  of  the  five  controllers.  The 
experiments  involve  position  trajectory  tracking  with  velocity  reversals  which 
exercise  the  problems  associated  with  friction  at  near-zero  velocities  and  the 
discontinuous  nature  of  friction  at  zero  velocity.  Tracking  performance  is  mea¬ 
sured  by  root-mean-square  (RMS)  position  error  which  provides  some  averaging 
of  external  effects.  Nonetheless,  unmodelled  effects  such  as  torque  ripple  and 
digital  sampling  rate  are  addressed  with  additional  experimentation. 

Section  2  of  this  paper  describes  the  different  models  of  friction  found  in  the 
literature  with  a  discussion  of  how  each  feature  of  friction  influences  servomech¬ 
anism  dynamics.  Section  3  discusses  friction-compensating  control  strategies 
and  presents  the  design  of  the  adaptive  and  estimation-based  controllers  tested 
experimentally.  Section  4  provides  the  details  of  the  experimental  program  and 
the  experimental  results.  Conclusions  along  with  some  suggestions  for  future 
work  are  given  in  Section  5. 


2  Friction  Structure  and  Dynamics 

Using  both  theory  and  experimentation,  researchers  in  a  number  of  fields  have 
developed  several  different  models  of  the  structure  and  dynamics  of  friction.  In 
selecting  a  friction  model  for  our  dynamic  friction-compensating  control  prob¬ 
lem,  it  is  important  to  consider  how  the  various  identified  features  of  friction 
influence  the  dynamics  of  a  servomechanism. 

Although  rolling  friction  is  a  physically  different  phenomenon  from  rubbing 
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Kins nc  plus 
Viscous  r fiction 


Figure  1:  (a)  Classical  Model  of  Static,  Kinetic  and  Viscous  Friction  and  (b) 
Stribeck  Friction 

or  sliding  friction,  the  models  discussed  below  attempt  to  describe  the  dynam¬ 
ics  of  a  system  with  rolling  or  sliding  friction.  Pure  rolling  friction  conditions 
occur  when  the  contact  between  two  surfaces  is  a  point.  However,  according  to 
[Rabinowicz,  1965],  the  contact  region  between  two  surfaces  is  typically  of  larger 
area  than  a  point  because  of  elastic  (and  possibly  plastic)  deformation  on  one  or 
both  of  the  surfaces.  The  resulting  “rolling”  friction  in  this  case  involves  a  com¬ 
bination  of  sliding  and  pure  rolling  friction.  In  fact,  although  the  sliding  velocity 
is  usually  small  compared  to  the  rolling  velocity,  sliding  friction  often  provides 
the  major  component  of  the  total  friction.  Consequently,  it  is  appropriate  to 
consider  the  same  models  for  sliding  friction  and  “rolling”  friction. 

Classical  friction  is  the  earliest  and  most  widely  used  model  of  friction.  The 
three  components  of  classical  friction:  kinetic  friction,  viscous  friction,  and 
static  friction,  are  illustrated  on  the  friction  versus  velocity  graph  of  Figure 
1(a).  Although  kinetic  friction  simply  provides  a  constant  retarding  force  to 
rubbing  surfaces,  it  also  introduces  a  discontinuity  at  zero  velocity.  As  a  result, 
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servomechanisms  performing  bi-directional  tasks  will  be  subject  to  the  discon¬ 
tinuity  during  every  velocity  reversal.  The  discontinuous  behavior  of  kinetic 
friction  can  be  classified  as  a  "'hard  nonlinearity'’ .  It  is  well-known  that  a  closed 
loop  system  with  a  hard  nonlinearity  can  produce  a  limit  cycle,  i.e. ,  self-sustained 
oscillations,  that  would  cause  poor  control  accuracy. 

Viscous  friction  results  from  the  viscous  behavior  of  a  fluid  lubricant  layer 
between  two  rubbing  surfaces.  As  shown  in  Figure  1(a),  viscous  friction  is 
represented  as  a  linear  function  of  velocity. 

Static  friction  is  the  force  required  to  initiate  motion  from  rest.  Typically, 
the  magnitude  of  static  friction  is  greater  than  the  magnitude  of  kinetic  friction 
which  can  lead  to  intermittent  motion  known  as  “stick-slip”.  Stick-slip  manifests 
itself  as  repeated  sequences  of  sticking  between  two  surfaces  with  static  friction 
followed  by  sliding  or  slipping  of  the  two  surfaces  with  kinetic  friction.  For  the 
servomechanism  control  problem,  stick-slip  can  diminish  control  accuracy.  The 
stick-slip  cycling  can  be  avoided  if  damping  and  stiffness  are  sufficiently  high. 

The  classical  lumped  friction  model  Ff  of  static  friction  Fs,  kinetic  friction 
Fk,  and  viscous  friction,  which  depend  on  the  applied  tangential  force  F,  velocity 
V ,  and  coefficient  of  viscous  friction  /i,  is  as  follows: 

_  |  Fksgn{V)  +  gV  if  V  ±  0  ^ 

{  Fssgn(F)  if  V  -  0 

Some  experimentalists  have  noted  that  in  machines  with  rubbing  parts  more 
complicated  and  numerous  than  a  single  body  sliding  over  a  second  body,  the 
magnitudes  of  kinetic,  viscous,  and  static  friction  are  not  the  same  in  the  pos¬ 
itive  and  negative  directions  [Canudas  et  al.,  1986,  Armstrong- Helouvry,  1991, 
Wang,  1987].  A  more  general  model  of  friction  that  accounts  for  this  asymmetry 
uses  different  friction  coefficients  in  the  positive  and  negative  directions. 
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Figure  2:  Pre-Sliding  Displacement  Phenomenon  (Dahl  Effect) 

Contrary  to  the  predictions  derived  from  the  classical  friction  model,  re¬ 
searchers  including  the  authors  of  [Courtney-Pratt  and  Eisner,  1957]  and  others 
have  found  that  small  relative  displacements  between  two  bodies  in  contact  do 
occur  when  the  applied  relative  tangential  force  is  less  than  the  static  friction. 
This  pre-sliding  displacement  phenomenon  is  illustrated  in  Figure  2  which  shows 
friction  Fj  as  a  function  of  displacement  x  based  on  experimental  results. 


In  [Dahl,  1977]  the  author  provided  a  model  of  the  pre-sliding  displacement 
phenomenon,  known  as  the  “Dahl  model”,  that  assumes  friction  Ff  is  a  function 
of  displacement  x  and  time  t  such  that 


dFf(x,t)  dF/(x,t)  .  dF/(x,t) 
—  •  x  + 


dtt 


dx 


dt 


with  dFf/dt  =  0,  and 


dFf(x,t)  .  Ft  • 


(3) 


<7  and  Ffc  are  as  shown  in  Figure  2,  and  i  is  an  exponent  parameter  that  Dahl 
empirically  derived  to  be  i  «  1.5. 


A  friction-compensating  adaptive  controller  based  on  Dahl's  model  was  de¬ 
signed  and  successfully  used  for  the  stabilization  of  an  airborne  pointing  and 
tracking  system  [Walrath,  1984],  Walrath  found  from  experimentation  that 
friction  responds  continuously  to  velocity  reversals.  Using  the  classical  discon¬ 
tinuous  static-kinetic  friction  model,  Walrath  could  not  re-create  this  smooth 
behavior.  Dahl’s  model,  on  the  other  hand,  predicted  the  expected  smooth 
behavior. 

The  classical  friction  model  predicts  behavior  more  characteristic  of  a  system 
that  spends  longer  periods  of  time  at  zero  velocity.  Specifically,  the  magnitude 
of  static  friction  is  dependent  on  the  length  of  time  the  surfaces  are  at  rest,  i.e. , 
the  “dwell  time”.  This  static  friction  dependence  on  dwell  time  explains  why 
under  stick-slip  conditions  the  amplitude  of  the  stick-slip  limit  cycle  is  observed 
to  decrease  with  increasing  velocity. 

While  the  simple  static  plus  kinetic  friction  model  offers  an  intuitive  expla¬ 
nation  for  the  possibility  of  stick-slip  oscillations,  it  inadequately  justifies  the 
existence  of  these  limit  cycles  for  all  the  conditions  under  which  they  have  been 
observed.  However,  several  researchers,  e.g.,  [Armstrong-Helouvry,  1991]  have 
experimentally  confirmed  that  friction  varies  with  velocity  as  depicted  in  Figure 
1(b).  The  implications  of  the  Stribeck  effect  with  regard  to  servomechanism 
dynamics  include  an  increased  likelihood  of  stick-slip  limit  cycling  at  low  ve¬ 
locities.  Of  the  many  empirical  models  derived  for  friction  incorporating  the 
Stribeck  effect,  the  following  is  the  most  popular: 

Ff(V)  =  Fksgn(V)  +  +  (Fs  -  Fk)e~^v^2 sgn(V)  .  (4) 

where  Vstr  is  the  critical  Stribeck  velocity. 

Frictional  lag  is  one  other  feature  of  friction  that  may  also  have  a  significant 
impact  on  dynamics.  While  we  do  not  consider  this  in  the  present  study,  we  note 
the  considerable  empirical  evidence  that  has  recently  become  available  indicating 


8 


that  friction  does  not  respond  instantaneously  to  a  change  in  velocity.  The  pri¬ 
mary  work  here  is  due  to  geophysicists  [Rice  and  Ruina,  1983]  who  use  stick-slip 
for  earthquake-related  predictions.  Hess  and  Soom  [Hess  and  Soom,  1990]  also 
found  strong  evidence  of  frictional  lag,  in  their  experiments  on  a  flat  steel  button 
rubbing  against  a  rotating  steel  disk.  Frictional  lag  makes  stick-slip  instabilities 
less  likely.  Because  a  decrease  in  friction  occurs  slowly  when  velocity  is  increased, 
stiff  systems  will  not  experience  stick-slip  [Armstrong-Helouvry,  1991].  We  hope 
to  investigate  these  aspects  in  future  work. 


3  Friction-Compensating  Control 

If  a  system  with  friction  is  linear  and  is  to  be  operated  only  at  relatively  high  ve¬ 
locities  without  changing  directions,  i.e.,  without  crossing  zero  velocity,  friction 
can  be  modelled  as  a  linear  function  of  velocity.  Under  these  conditions,  stan¬ 
dard  PID  design  techniques  can  be  applied  to  the  dynamics  of  the  linear  system 
plus  viscous  friction  with  reliable  results.  On  the  other  hand,  if  the  system  is 
to  be  operated  at  low  velocities  or  with  direction  reversals,  then  the  standard 
PID  design  techniques  may  be  unsuitable  and  tracking  accuracy  may  prove  in¬ 
adequate.  Additionally,  to  prevent  limit  cycling  due  to  the  static-kinetic  friction 
discontinuity  at  zero  velocity  or  the  Stribeck  effect,  a  PID  controller  must  have 
sufficiently  high  “damping”  Kd  and  “stiffness”  Kp.  However,  high  gain  control 
has  its  own  practical  disadvantages  such  as  introducing  instability  in  a  compliant 
drive  train  or  saturating  an  actuator. 

One  popular  alternative  is  the  use  of  a  technique  such  as  dither  which  trans¬ 
forms  the  dynamics  of  a  system  with  a  discontinuity  into  smooth  dynamics  that 
can  be  more  easily  controlled  with  standard  techniques.  Dither  is  a  high  fre¬ 
quency  signal  added  to  the  error  signal  in  a  feedback  loop  before  it  is  input  to 
the  system.  If  the  frequency  is  chosen  to  be  higher  than  the  cut-off  frequency 
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of  the  system,  the  high-frequency  behavior  is  filtered  out  leaving  only  the  low- 
frequency  "average”  response. 

Pulse-width  modulation  is  another  commonly  used  and  effective  smoothing 
technique  that  also  works  on  the  principle  of  averaging.  However,  both  dither 
and  pulse-width  modulation  have  inherent  disadvantages.  For  example,  analy¬ 
sis  and  prediction  of  system  characteristics  such  as  stability  and  robustness  are 
difficult  to  perform  when  dither  or  pulse-width  modulation  is  applied.  Addi¬ 
tionally,  dither  can  cause  mechanical  problems  in  a  system  such  as  a  robot  by 
exciting  vibrations. 

Friction-compensating  adaptive  and  estimation-based  controllers  are  nonlin¬ 
ear  controllers  that  can  be  designed  to  take  advantage  of  what  is  known  about 
the  structure  of  friction.  Adaptive  control  strategies,  in  particular,  are  nat¬ 
urally  suited  to  the  problem  of  friction  compensation  because  they  generate  a 
time- varying  control  law  that  tracks  slowly  varying  system  parameters,  and  they 
provide  system  identification  when  an  accurate  system  model  is  not  available. 
Three  adaptive  and  estimation-based  controllers  are  investigated  in  this  paper. 

For  reference,  the  dynamics  of  the  dc  motor  used  in  this  paper  are  described 

by 

0p(t)  +  ci$p(t)  =  -c2Tj  +  c3u(t)  ,  (5) 

where  0P,  9p:  9P  are  plant  angular  position,  velocity,  and  acceleration,  Tf  is  the 
friction  and  may  depend  on  9P,  9P,  etc.,  u  is  the  control  input,  and  ci,  c2,  C3  are 
constants  (ci  includes  viscous  friction). 

The  first  of  the  adaptive  and  estimation-based  controllers  is  an  adaptive 
controller  (referred  to  as  AEC  I)  based  on  the  Model  Reference  Adaptive  Con¬ 
trol  (MRAC)  approach  proposed  by  [Gilbart  and  Winston,  1974]  in  their  work 
on  control  of  a  satellite-tracking  telescope.  However,  while  Gilbart  and  Win¬ 
ston  considered  velocity  trajectory  tracking  in  their  design,  the  controller  in  the 
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present  paper  was  required  to  handle  position  trajectory  tracking.  As  a  result, 
AEC  I  is  a  modified  version  of  the  Gilbart  and  Winston  design,  incorporating 
filters  to  reduce  the  order  of  the  system. 


Additionally,  Gilbart  and  Winston  assumed  the  classical,  symmetrical  kinetic 
plus  viscous  friction  model  in  the  dynamic  system  equations.  To  accommodate 
the  observed  asymmetrical  nature  of  friction  in  the  motor  [Wang,  1987],  kinetic 
friction  is  modelled  in  AEC  I  as  follows: 


Tf  =  oq 


'  sgn(6p)  +  l' 


+  0-1 


sgn(6p)  —  1 


(6) 


where  oq  and  a 2  represent  the  magnitude  of  kinetic  friction  in  the  positive  and 
negative  directions,  respectively. 


While  Gilbart  and  Winston  used  only  proportional  feedback  control  in  ad¬ 
dition  to  the  adaptive  control,  our  controller  uses  a  control  input  based  on 
the  computed  torque  method  with  integration  in  addition  to  the  adaptive  con¬ 
trol.  Figure  3  shows  the  block  diagram  for  AEC  I.  A  stability  analysis  of  this 
controller  which  follows  [Gilbart  and  Winston,  1974]  and  employs  Lyapunov’s 
direct  method  yields  the  result  that  (e,  e)  =  (0, 0)  is  an  asymptotically  sta¬ 
ble  equilibrium  point  where  e  =  9m  —  9P,  e  =  9m  —  9P  and  9m  and  9m  are 
the  ideal  model  velocity  and  position,  respectively.  Details  are  to  be  found  in 
[Ehrich,  1991]. 

The  second  controller  (referred  to  as  AEC  II)  is  based  on  the  method  devel¬ 
oped  by  [Craig,  1988]  in  his  design  of  an  adaptive  robotic  manipulator  controller. 
As  in  the  work  of  Gilbart  and  Winston,  Craig  used  an  MRAC  approach  and  as¬ 
sumed  the  classical,  symmetrical  model  of  kinetic  plus  viscous  friction.  Four 
different  versions  of  the  AEC  II  controller  are  developed,  each  with  a  different 
model  of  friction.  For  reference,  we  solve  equation  (5)  for  input  u  and  rewrite 
as 

u  =  (1  /c3)0p  +  (c!/c3)0p  +  (c2/c3)T/  =  ( 1  / c3 ) 9p  +  Q  .  (7) 
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Figure  3:  AEC  I 


Friction  model  (a)  is  the  simple  kinetic  plus  viscous  friction  model  used  by 
Craig  in  his  design  such  that 


Q{a)  =  pJp  +  P2 sgn(9p)  , 


(3) 


where  p\  and  p2  are  the  unknown  parameters.  Model  (b)  assumes  the  asymmet¬ 
ric  kinetic  plus  viscous  friction  model  represented  by 

Q(b)  =  Plopi±is^m  +  P2epi^M 

i  r,  sgn(ep)+l  ,  sgn(flp)-l 
'  2  ‘  ^4  2  ’ 

where  pi  and  p2  are  the  unknown  viscous  friction  parameters  in  the  positive  and 
negative  directions,  respectively  and  p3  and  p4  are  the  unknown  kinetic  friction 
parameters  in  the  positive  and  negative  directions,  respectively. 


Model  (c)  includes  a  linear  model  of  Stribeck  friction  in  addition  to  kinetic 


plus  viscous  friction.  Stribeck  friction  can  be  modelled  according  to  (4)  for 


12 


Figure  4:  AEC  II 

V  =  9P  and  Vstr  —  OstT.  However,  AEC  II  requires  dynamics  that  are  linear  in  the 
unknown  parameters.  To  derive  a  linear-in-parameters  model,  the  exponential 
in  (4)  is  replaced  by  its  Taylor  series  approximation  which  yields 

Q(c)  =  piK  +  p2sgn(dp )  +  p30p2 sgn(6p)  .  (10) 

For  comparative  purposes,  Model  (d)  includes  kinetic  plus  vis¬ 
cous  plus  Stribeck  friction  terms  according  to  the  linearized  model  of 
[Canudas  de  Wit,  1989]  such  that 

Q{d)  =  Pidp  +  P2 sgn(9p)  +  P3\9p\1/2sgn(6p)  .  (11) 

The  block  diagram  for  AEC  II  is  shown  in  Figure  4.  A  stability  analysis  that 
follows  [Craig,  1988]  and  uses  the  Kalman- Yakubovich  Lemma  and  Lyapunov’s 
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direct  method  shows  that  (e,  e)  =  (0,0)  where  e  =  O.j  —  6P  and  e  =  9d  —  0P  (see 
[Ehrich,  1991]). 


The  third  controller,  AEC  III,  is  an  estimation-based  controller  that  follows 
the  work  of  [Walrath,  1984]  described  in  the  previous  section.  From  experimen¬ 
tal  results,  Walrath  postulated  the  following  first-order  model  of  the  bearing 
friction  Tj: 

r^-  +  T,  =  Tcsgn(0r)  (12) 

where  Tc  is  the  rolling  bearing  friction  (i.e.  kinetic  friction)  and  r  is  a  time 
constant.  Clearly  (12)  is  a  rewriting  of  the  Dahl  friction  model  of  (3)  with  i  —  1 
and  replacing  Fj  with  T/,  Fjc  with  Tc  and  x  with  dp  where 


r  = 


Tc 

a\dp\ 


(13) 


Walrath  incorporated  the  friction  model  (12)  into  his  controller  by  using 
it  to  predict  friction  torque.  He  empirically  determined  the  value  of  r  for  a 
given  control  experiment  by  repeating  the  experiment  many  times,  each  time 
varying  only  r,  until  the  optimum  r  was  found,  i.e.,  the  r  that  yielded  the 
minimum  stabilization  error.  Based  on  optimum  r  values  (ropt)  calculated  for 
a  range  of  operating  conditions,  Walrath  empirically  derived  a  linear  prediction 
of  ropi  as  a  function  of  the  inverse  of  the  RMS  system  acceleration  9rms.  This  is 
inconsistent  with  Dahl’s  model  which  predicts  that  r  is  inversely  proportional 
to  velocity  (13). 


On  the  other  hand,  a  consistent  relationship  for  r  was  found  for  the  electric 
motor  of  the  present  study.  Experiments  similar  to  Walrath’s  were  performed 
on  the  motor  to  determine  r  as  a  function  of  operating  conditions.  Figure  5(a) 
shows  the  results  of  one  of  the  experiments.  RMS  angular  position  error 
is  plotted  against  u  =  1/r  for  sinusoidal  trajectory  tracking  with  trajectory 
frequency  /  —  0.25  Hz  and  amplitude  A  =  0.25  rad.  Each  point  on  the  plot 
represents  a  single  repetition  (320  seconds  long)  of  the  experiment.  The  optimal 
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Figure  5:  (a)  Experimental  Results  Used  to  Determine  Optimal  u  (/  —  0.25 
Hz)  and  (b)  Optimal  Values  of  cu  as  Function  of  Operating  Conditions 


u  ( u)opt )  was  selected  to  correspond  to  the  minimum  RMS  error.  Figure  5(b) 
plots  uopt  for  all  the  experiments  as  a  function  of  RMS  velocity  and  shows  uopt 
to  be  a  linear  function  of  RMS  velocity. 


Additionally,  while  Walrath’s  controller  only  used  proportional  feedback  con¬ 
trol  in  conjunction  with  the  friction  compensation,  AEC  III  uses  feedback  control 
based  on  the  computed  torque  method  with  integration.  Figure  6  shows  a  block 
diagram  of  AEC  III. 

Since  Walrath  did  not  provide  a  stability  analysis  for  his  controller,  an  orig¬ 
inal  stability  analysis  has  been  developed  based  on  the  passivity  formalism. 
Consider  AEC  III  with  derivative  feedback  only.  Let  the  forward  part  of  the 
loop  including  the  motor,  friction  prediction  and  derivative  gain  K,i  be  repre¬ 
sented  by  H\.  The  feedback  part  of  the  loop  is  represented  by  H2  which  implies 
that  H2  =  1  from  Figure  6.  Then  it  follows  from  [Hill  and  Moylan,  1977],  that 
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Figure  6:  AEG  III 

if  both  H\  and  H2  are  Y-strongly  passive  (YSP),  the  feedback  system  is  asymp¬ 
totically  stable.  The  fact  that  H2  is  SPR  implies  that  H2  is  YSP.  It  remains  to 
prove  that  Hi  is  YSP. 


System  H\  can  be  described  by  the  following  equations: 

9p  +  c\9p  —  —c2Tj  +  c2Tf  +  c3Kde  ,  (Id) 

Tf  =  -b\6p\Tj  +  bTjp ,  (15) 

tf  =  -b\dp\tf  +  bTcQp.  (16) 

where  the  A  indicates  a  predicted  value  and  we  have  assumed  u  —  b\6p\  and 

Cj  =  b\9p\. 


Define  the  state  (aq,  x2,  x3,  x4)  =  (9p,9p,Tf,Tf)  and  let  u  —  e.  Now  define 
(f>  i 1Z4  i — >■  72.  cts 


4>{x) 


c3kiq[^x +  X3(xi  -  7i/&)  —  \Tcb(x i  -  71  /b)2  +  f  bxix3\x2\ 
+x4(xi  —  72 /b)  —  \Tcb(x i  —  72/S)2  +  f  bx\Xi  j x2 1] 

(17) 
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where  K,  71}  and  72  are  constants  to  be  defined  later.  Differentiating  (IT)  with 
respect  to  time  yields  the  supply  rate 

cf>(x)  —  uy  —  ey2  .  (18) 

where  e  =  Cijc^Ks.  >  0  and  7 1?  72,  and  K  are  chosen  appropriately  for  the  two 
cases  X2  =  dp  >  0  and  x2  =  0P  <  0. 

For  x2  —  9P  >  0  (or  x2  =  9V  <  0),  0:10:3  =  Tf9p  >  0  and  =  T/#p  > 
0.  Additionally,  07  and  0:3  are  both  bounded  since  0  <  x\  —  9V  <  27r  and 
— Tc  <  £3  =  Tf  <  Tc.  Therefore  d>(-)  from  (17)  is  bounded  below.  Define  C  as 
the  greatest  lower  bound  of  d>(x).  Then  there  exists  some  x0  6  7T.4  such  that 
<fi(xo)  =  C.  Next  define  <p  :  1Z4  1Z  such  that 

4>{y)  =  <p(y  +  *0)  -  c  .  (19) 

Then  <^(0)  =  0  and  <j>(y)  >  0  Vy  and  <f>(y)  =  <j>(y).  Thus,  according  to 
[Hill  and  Moylan,  1977]  H\  is  YSP  and  the  feedback  system  of  AEC  III  is  input- 
output  asymptotically  stable.  This  implies  that  9P  will  follow  9 


4  Experimental  Program  and  Results 

The  experimental  system  consisted  of  a  direct-drive,  brush-type  dc  motor,  an¬ 
gular  position  and  velocity  sensors,  a  power  amplifier,  an  IBM  AT  PC,  and 
supporting  hardware  and  software  for  communication  and  control.  The  IBM 
PC  was  used  to  control  the  operation  of  the  motor.  A  50  Hz  sampling  rate  was 
used  throughout  the  experimental  program. 

Angular  position  of  the  motor  was  measured  by  a  12  bit  absolute  optical 
shaft  encoder.  The  position  data  was  transmitted  to  the  PC  with  a  measure¬ 
ment  resolution  of  0.00154  radians  (0.088  degrees).  A  tachometer  measured  the 
angular  velocity  of  the  motor  with  a  resolution  of  0.012  rad/s. 
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The  motor  system  was  modelled  according  to  (5)  where  u  is  the  input  voltage 
to  the  motor.  For  system  simulations,  friction  was  modelled  as  symmetric  kinetic 
plus  viscous  friction.  The  friction  parameters,  assumed  to  be  constant,  were 
measured  in  previous  work  by  [Wang,  1987].  The  motor  model  was  successfully 
verified  by  comparing  results  of  voltage  pulse  experiments  on  the  motor  to  results 
of  identical  simulated  experiments  on  the  motor  model. 

The  motor  model  along  with  a  simulation  of  a  standard  PID  controller  was 
then  used  to  design  an  optimized  benchmark  PID  controller  for  position  tracking. 
The  PID  controller  gains  were  selected  to  optimize  the  motor  response  to  a  0.25 
radian  step  demand  in  angular  position.  CONSOLE,  a  numerical  optimization 
tool  described  in  [Fan  et  ah,  1987],  was  used  to  perform  the  optimization.  Two 
functional  objectives  were  specified,  one  to  minimize  the  overshoot  of  the  step 
response  and  one  to  maximize  the  rise  of  the  step  response.  The  effectiveness  of 
the  optimized  gains  was  verified  by  a  0.25  radian  step  experiment  on  the  motor. 

The  experimental  program  on  the  electric  motor  consisted  of  comparative 
position  trajectory  tracking  tests  using  five  different  controllers:  AEC  I,  II, 
III,  a  controller  with  dither  and  the  optimized  benchmark  PID  controller.  The 
controller  with  dither  was  implemented  identically  to  the  PID  controller  except 
that  a  dither  signal  was  added  to  the  control  input.  The  frequency  of  the  dither 
signal  was  25  Hz  which  is  the  maximum  possible  given  the  50  Hz  sampling  rate. 
The  amplitude  of  the  dither  signal  was  more  than  twice  the  magnitude  of  the 
static  friction. 

In  each  experiment  the  motor  was  required  to  track  a  sinusoidal  position 
trajectory  such  that: 

9d  =  Asin(2ir  ft)  (20) 

where  A  is  the  amplitude  and  /  the  frequency  of  the  demanded  trajectory. 
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This  required  sinusoidal  motion  provided  a  useful  means  for  investigating  fric¬ 
tion  compensation  since  the  motor  was  forced  to  repeatedly  pass  through  zero 
velocity  where  friction  behavior  is  most  difficult  to  control.  The  sinusoidal  mo¬ 
tion  also  provided  a  reasonably  realistic  scenario  since  manipulators  are  often 
required  to  perform  repetitive  tasks  that  demand  sinusoidal  joint  motions. 

The  sinusoidal  trajectories  tracked  in  the  experimental  program  ranged  in 
frequency  /  from  0.1  Hz  to  1.0  Hz.  The  lower  limit  of  this  range  was  selected  to 
minimize  motor  velocity  and  to  avoid  large  errors  due  to  velocity  measurement 
resolution.  The  upper  limit  of  this  range  was  selected  to  maximize  motor  velocity 
without  generating  gross  errors  due  to  the  limitations  of  the  50  Hz  sampling 
rate.  A  0.25  radian  amplitude  A  was  used  for  all  sinusoidal  trajectory  tracking 
experiments. 

The  RMS  position  error  for  each  of  the  controllers  and  each  experiment  is 
listed  in  Table  1.  Each  RMS  position  error  is  calculated  based  on  16  seconds  (800 
samples)  of  data.  As  indicated  in  Table  1,  the  controller  with  dither  does  not 
significantly  improve  the  tracking  performance  as  compared  to  the  benchmark 
PID  controller.  This  is  due  to  the  fact  that  the  dither  frequency  is  limited  to 
25  Hz.  AEC  I,  II,  and  III,  on  the  other  hand,  all  effectively  improved  tracking 
performance  for  the  range  of  sinusoidal  trajectory  frequencies  tested.  The  lower 
values  of  percent  reduction  in  RMS  error  for  the  0.1  Hz  experiment  are  most 
likely  due  to  the  resolution  of  the  position  measurements. 

Figures  7,  8,  and  9  show  the  absolute  value  of  PID  controller  position  error 
compared  to  the  position  error  from  AEC  I,  II,  and  III,  respectively.  In  each 
case  the  PID  controller  result  is  shown  in  a  solid  line  and  the  nonlinear  controller 
result  is  shown  in  a  dashed  line.  Note  that  the  two  seconds  (100  samples)  shown 
correspond  to  the  first  two  seconds  of  the  experimental  results  beginning  at  the 
time  listed  in  the  fourth  column  of  Table  1. 
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Type 

Trajectory 

RMS  Position 

Time  Elapsed 

%  Error 

of 

Frequency  / 

Error 

Before  Error 

Reduction 

Controller 

(Hz) 

(rad) 

Calculated  (s) 

from  PID 

PID 

1.0 

0.0106 

0 

- 

Dither 

1.0 

- 

- 

- 

AEC  I 

1.0 

0.0048 

400 

55 

AEC  II 

1.0 

0.0079 

12 

25 

AC  III 

1.0 

0.0054 

12 

49 

PID 

0.5 

0.0069 

0 

- 

Dither 

0.5 

0.0066 

0 

4 

AEC  I 

0.5 

0.0033 

400 

52 

AEC  II 

0.5 

0.0041 

12 

41 

AEC  III 

0.5 

0.0028 

12 

59 

PID 

0.25 

0.0063 

0 

- 

Dither 

0.25 

0.0055 

0 

13 

AEC  I 

0.25 

0.0043 

400 

32 

AEC  II 

0.25 

0.0037 

12 

41 

AEC  III 

0.25 

0.0036 

12 

43 

PID 

0.1 

0.0060 

- 

Dither 

0.1 

0.0045 

25 

AEC  I 

0.1 

- 

- 

- 

AEC  II 

0.1 

0.0044 

12 

27 

AEC  III 

0.1 

0.0040 

12 

33 

Table  1:  Results  of  Sinusoidal  Tracking  Experiments 

A  comparison  of  these  plots  and  the  performance  results  of  Table  1  suggests 
that  AEC  II  is  not  as  effective  at  friction  compensation  as  AEC  III.  This  can 
be  explained  by  noting  that  there  are  a  couple  of  relatively  large  error  peaks 
in  Figure  9  for  AEC  II.  These  occur  because  AEC  II  overcompensates  when 
friction  changes  instantaneously.  The  large  error  peaks  seen  in  Figure  11  cor¬ 
respond  in  time  to  the  instantaneous  friction  changes.  This  overcompensation 
and  corresponding  large  error  indicates  that  the  classical  friction  model  does  not 
describe  friction  during  transient  velocity  reversals  as  well  as  the  Dahl  model. 

The  numbers  listed  in  Table  1  are  averages  of  results  from  experiments  made 
over  a  period  of  a  few  days.  However,  over  the  course  of  about  six  months 
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.Figure  7:  Position  Error  for  AEG  I  Compared  to  Benchmark  PID  Controller 

during  which  these  experiments  were  performed,  there  was  a  great  deal  of  re¬ 
peatability  in  the  percent  reduction  in  RMS  position  error  achieved  by  the  two 
adaptive  controllers.  The  experiments  were  run  during  different  seasons  and 
during  different  stages  of  motor  “warm-up”  such  that  friction  parameters  may 
have  varied  from  experiment  to  experiment  due  to  temperature  differences.  Ad¬ 
ditionally,  over  the  six  month  period  the  friction  parameters  may  have  changed 
due  to  system  aging.  The  fact  that  the  adaptive  controllers  were  consistently 
effective  under  these  varying  conditions  provides  evidence  for  the  effectiveness 
of  the  adaptability  of  these  controllers. 

The  results  for  AEC  II  provided  in  Table  1  correspond  to  experiments  per¬ 
formed  using  Model  (a)  friction.  Experiments  using  Model  (b)  friction  showed 
an  insignificant  change  in  performance.  This  result  indicates  that  for  this  set 
of  experimental  conditions,  additional  adaptive  terms  to  account  for  friction 
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Figure  8:  Position  Error  for  AEC  II  Compared  to  Benchmark  PID  Controller 

asymmetries  are  not  necessary.  Experiments  using  Models  (c)  and  (d)  for  fric¬ 
tion  showed  only  a  slight  improvement  in  performance  as  compared  to  the  Model 
(a)  experiments.  Since  the  additional  friction  terms  in  Model  (c)  and  Model  (d) 
were  intended  to  account  for  Stribeck  friction,  this  result  can  be  attributed  to 
the  fact  that  Stribeck  friction  is  probably  not  completely  measured  by  the  ex¬ 
perimental  system  since  the  critical  Stribeck  velocity  may  be  approximately  of 
the  same  or  lower  order  as  the  velocity  measurement  resolution. 

The  main  disadvantage  associated  with  AEC  I  was  its  excessively  slow  rate 
of  adaptation.  Although  tracking  performance  began  to  improve  immediately, 
it  was  not  until  400  seconds  into  the  experiment  that  the  best  results  were 
achieved. 

AEC  II  had  the  disadvantage  that  for  the  implementation  that  yielded  the 
best  results,  the  adaptive  parameters  tended  to  drift  and  performance  deteri¬ 
orated  after  a  while.  This  could  be  avoided  by  resetting  the  parameters  when 
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Figure  9:  Position  Error  for  AEC  III  Compared  to  Benchmark  PID  Controller 
they  went  out  of  a  predefined  range  as  suggested  by  [Craig,  1988]. 

AEC  III,  on  the  other  hand,  was  very  reliable  and  performed  best  of  all 
the  controllers.  However,  this  controller  is  at  a  disadvantage  in  that  it  requires 
a  lengthy  experiment  up  front  to  determine  the  constants  in  the  relationship 
between  the  friction  time  constant  r  and  the  RMS  velocity.  Additionally,  since 
kinetic  friction  Tc  is  held  constant,  AEC  III  is  not  best  suited  for  adapting  to 
changes  in  friction  due  to  temperature  or  aging.  This  could  be  fixed,  however, 
by  adding  an  adaptive  component  to  update  Tc. 

Finally,  torque  ripple  in  the  motor  adds  a  position-dependent  component 
to  the  motor  dynamics  and  could  have  affected  how  each  of  the  controllers 
performed.  All  of  the  data  in  Table  1  applies  to  experiments  run  such  that  the 
initial  position  was  0.0  radians.  However,  it  was  observed  that  the  performance 
of  the  controllers  varied  when  different  initial  positions  were  used.  To  investigate 
this  torque  ripple  effect,  the  experiment  with  /  =  0.5  Hz  was  run  again  on 
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Initial 

RMS  Position  Error 

%  RMS  Position  Error 

Position 

with  PID  Controller 

Reduction  from  PID 

(rad) 

(rad) 

AEC  11 

AEC  III 

-2.698 

0.0086 

16 

23 

-2.484 

0.0082 

26 

33 

-2.075 

0.0091 

21 

23 

-1.546 

0.0081 

30 

35 

-1.080 

0.0079 

30 

39 

-0.709 

0.0079 

28 

29 

-0.261 

0.0072 

33 

36 

-0.069 

0.0071 

32 

41 

0.000 

0.0067 

34 

52 

0.086 

0.0068 

34 

46 

0.689 

0.0060 

28 

42 

1.172 

0.0061 

15 

56 

1.758 

0.0066 

IS 

39 

2.484 

0.0076 

32 

37 

2.720 

0.0074 

22 

39 

Average 

0.0074 

27 

38 

Std.  Dev. 

0.00087 

6.4 

CO 

Table  2:  Experimental  Results  at  Different  Initial  Positions  (/  =  0.5  Hz) 

the  PID  controller  and  AEC  II  and  III  at  15  different  initial  positions  chosen 
randomly.  The  results  of  these  experiments  are  shown  in  Table  2.  According 
to  these  results,  neither  AEC  II  nor  III  performed  on  average  as  well  as  at  an 
initial  position  of  0.0  radians.  However,  the  greater  effectiveness  of  AEC  III 
relative  to  AEC  II  was  observed  at  every  initial  position.  From  Table  2,  it  can 
be  concluded  that  while  torque  ripple  does  affect  somewhat  the  performance  of 
these  two  adaptive  controllers,  it  does  not  affect  their  relative  effectiveness. 

The  effect  of  digital  sampling  rate  on  the  performance  of  AEC  II  and  III  was 
also  investigated  by  repeating  the  experiments  of  the  experimental  program  with 
a  100  Hz  sampling  rate.  Table  3  lists  the  results  of  these  experiments.  According 
to  Table  3,  the  increased  sampling  rate  did  not  have  a  dramatic  effect  on  the 
performance  of  AEC  III.  However,  AEC  II  performed  significantly  better  with 
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Trajectory 
frequency  / 
(Hz) 

%  RMS  Position  Error  Reduction  from  PID 

AEC  II 

AEC  III 

50  Hz 

100  Hz 

50  Hz 

100  Hz 

1.0 

25 

64 

49 

42 

0.5 

41 

60 

59 

56 

0.25 

41 

45 

43 

56 

0.1 

27 

35 

33 

29 

Table  3:  Comparison  of  Tracking  Experiment  Results  with  50  Hz  and  100  Hz 
Sampling  Rates 

the  100  Hz  sampling  rate  than  with  the  50  Hz  sampling  rate,  particularly  for 
experiments  with  /  =  0.5  Hz  and  /  =  1.0  Hz.  This  improved  performance 
may  be  explained  by  the  fact  that  overcompensation  provided  by  AEC  II  for 
instantaneous  changes  in  friction  is  not  as  prolonged  with  a  100  Hz  sampling 
rate  as  it  is  with  a  50  Hz  sampling  rate.  Based  on  the  results  of  Table  3,  one 
can  conclude  that  the  relative  effectiveness  of  AEC  II  and  III  is  greatly  affected 
by  the  digital  sampling  rate. 


5  Conclusions 

In  this  paper  a  comprehensive  investigation  was  presented  of  control  strategies 
for  friction  compensation  in  servomechanisms  performing  low-velocity  position 
tracking.  The  major  conclusions  of  the  investigation  are  as  follows: 

•  AEC  I,  II,  and  III  all  provide  improved  servomechanism  position  con¬ 
trol  compared  to  an  optimized  PID  controller  and  controller  with  (limited 
frequency)  dither  for  low-frequency  sinusoidal  position  trajectory  tracking 
experiments  on  a  direct-drive  dc  motor.  Additionally,  the  experimental  re¬ 
sults  of  this  paper  coupled  with  the  results  of  [Gilbart  and  Winston,  1974, 
Craig,  1988,  Walrath,  1984]  provide  evidence  for  the  general  applicability 
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of  these  adaptive  and  estimation-based  controllers. 


•  The  Dahl  model  provides  a  realistic  and  reliable  model  of  friction,  partic¬ 
ularly  during  sinusoidal  motion  of  the  mechanism.  Evidence  for  this  can 
be  found  (1)  in  the  fact  that  the  empirically  derived  model  of  the  friction 
time  constant  r  as  a  linear  function  of  velocity  is  consistent  with  Dahl’s 
original  model  and  (2)  by  the  relatively  high  effectiveness  of  AEC  III  which 
is  based  on  the  Dahl  model.  This  conclusion  is  noteworthy  since  friction 
is  typically  considered  to  behave  according  to  the  classical  friction  model. 

•  Mechanical  considerations  such  as  torque  ripple  and  digital  sampling  rate 
play  an  important  role  in  the  performance  of  the  adaptive  and  estimation- 
based  controllers. 

Further  research  should  be  pursued  to  understand  the  relationship  between 
the  classical  friction  model  and  the  Dahl  friction  model.  A  determination  of  how 
to  link  the  Dahl  model  of  pre-sliding  displacements  with  the  classical  model  of 
sticking  and  sliding  would  provide  a  more  complete  and  cohesive  understand¬ 
ing  of  friction  that  could  potentially  be  used  to  improve  friction-compensating 
control  strategies. 
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